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Progress 


Atmospheric eradient studies: 


Work is continuing on the study of atmospheric gradients. We include a preprint 
outlining some recent results which have been submitted to the journal 
"Geophysical Research Letters". 


Solid earth tides: 


Work has begun on a study of solid Earth tidal deformations using the VLBI data 
set. At present we have examined deformations at the semi-diurnal tidal period 
using the IRIS data set. We presented preliminary results at the spring 1992 
meeting of the American Geophysical Union in Montreal, Canada. 
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Abstract . Spatial and temporal fluctuations in the refractive index of air will ef- 
fect estimates of atmospheric parameters determined from data obtained from radio 
signals which have propagated through the atmosphere. I present an expression for 
the covariance of a general vector of atmospheric parameters, and use this expres- 
sion to calculate the effects on estimated gradient parameters of a fro2en Kolmogorov 
turbulence field moving over a site at constant velocity. Numerical calculations are per- 
formed to investigate the implications for three techniques: ground-based microwave 
radiometry, and geodesy with the Global Positioning System and very long baseline 
interferometry. The results indicate that care must be taken in comparing gradient 
parameters determined using data from these different techniques. 

1. Introduction 

Atmospheric water vapor, located primarily in the lower troposphere, significantly 
effects the propagation of radio waves in this region by retarding and refracting them. 
The additive contribution of water vapor to the atmospheric refractive index is the wet 
refractivity. The wet refractivity varies spatially in an apparently random manner, 
with significant fluctuations occurring on a range of spatial scales, from microscopic 
scales (1 mm or less) to scales of several thousand kilometers. These fluctuations 
effect estimates of atmospheric parameters determined from observations sensitive to 
the atmospheric refractive index, such as ground-based remote microwave sensing of 
the atmosphere and space geodetic measurements. In particular, we will examine 
the effect of these fluctuations on the estimation of atmospheric-gradient parameters, 
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which have been investigated using both of these techniques [Davis et al, 1992; Rogers 
et al., 1992], 

The development of a statistical model for gradient parameters will follow closely 
the methods and formalism of Treuhaft and Lanyi [19S/] (hereafter referred to as 
TL). In particular, we will assume that the spatial fluctuations in the wet refractivity 
discussed above are described by Kolmogorov turbulence theory, and that temporal 
wet refractivity fluctuations are caused by the “frozen field” of spatial fluctuations 
borne horizontally by the wind. The assumption of Kolmogorov turbulence leads to 
the 2/3 law” for the structure function of wet refractivity fluctuations [Tatarskii, 
1961]: 

D X (R) = <(*(r + R) - x (r)) 2 ) = C 2 R 2 / 3 (1) 

v here r is a three-dimensional location vector, R a displacement vector (and R its 
length), x the wet refractivity, and D x the structure function. The parameter C 
governs the power of the fluctuations. Like TL, we will assume that C does not vary 
spatially or temporally. 

2. A Statistical Model for Estimated Parameters 

In quantifying the effect of refractivity fluctuations on estimated gradient pa- 
rameters, we must explicitly consider the method used to obtain the estimates. The 
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techniques mentioned above are not sensitive to fluctuations in the refractive index di- 
rectly, but instead are sensitive to the wet delay, which is the integral of the refractivity 
along the line of sight: 

rh esc € 

T (,6,t — / ds x(r(e,9,z),t) (2) 

Jo 

where r is the delay along the line-of-sight at time t in the direction of azimuth angle 
8 and elevation angle e, and ds is the differential of distance along the line-of-sight 
path. In (2), we have parametrized the position vector r in terms of the elevation and 
azimuth angles, and the vertical coordinate r. The integral extends in height to h, the 
effective height of the troposphere, after the method of TL. 

For an isotropic troposphere, the wet delay in (2) varies with elevation angle 
approximately as the cosecant of elevation angle. In order to compare propagation 
delays at different elevation angles, it is customary to use the “equivalent” zenith 
delay r : given by 

Ke,t = sin e = [ dz *(r(e, 8, z), t ) (3) 

Jo 

In (3), even though the integration is performed in the zenith direction, the location 
vector r still varies along the true line of sight. 

For ground-based microwave radiometry, the equivalent zenith delay forms our 
observable for the estimation of gradient parameters [Davis et al., 1992]. The entire 
set of observables can be represented by a vector y, where 


y. = 




( 4 ) 
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and the estimated parameters p, which include gradient parameters, can be estimated 
using the equation 

p = (A T A)-'A T y (5 ) 

where A is the so-called design or partial derivative matrix. The form (5) is simply 
the least-squares estimate of p assuming equal weight (i.e., equal variances) for the 
observations y. Although in general we may not have equal weights, the assumption 
simplifies calculations without changing the conclusions of this study. Furthermore, 
we will assume for purposes of this study that the observations y are free from ob- 
servational errors. This last assumption notwithstanding, the vector p is a random 
vector because the observations y are stochastic, as discussed above. Although Kol- 
mogorov turbulence is in general nonstation ary, during the period of observations we 
can define the average wet delay (y), and the deviation from the average Sy = y — (y) . 

The estimates of the parameters will then deviate from their “average” values by the 
amount 

<5p = {A T A)~ l A T 6y (6) 

The expression (6) can be used to determine statistics for Sp. The expected value of 
this vector is just 

(<5p) = (A T A) -1 A T {6y) = 0 ( 7 ) 

Its covariance matrix is therefore 

<Mp T ) = (A T A)- 1 A T {8ySy T )A(A T A)~ l (S) 
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where we have used the fact that (-4^.4) 1 is symmetric. Denoting the covariance 
matrix of 8y in (8) as T, the element T,y is, using (2)-(4), given by 


T ij = (* r «* A,*.- * T c* = / dz [ dz> ( S x(r(e i ,e i ,z),t i )6x(r(e j ,e j ,z'),t j )) (9) 

Jo Jo 

where 6 again denotes deviation from the average. The expectation in (9) is difficult 
to evaluate, because it involves refractivity separated not only in time but in space. 
However, using the frozen flow hypothesis discussed above, we have 




*')» t j) = $x(r (<b, *') ~ vAty, f.) 


( 10 ) 


where At,y = tj — ti and v is the (horizontal) wind velocity vector, which we will take 
to be constant. Using (10), (9) can be written 


T H= f dz [ dz> (tix(r(e i ,8 i ,z),t i )6x(r(€ j ,6 j ,z’)-vAt ij ,t i )) 
Jo Jo 

= f dz [ dz' K x (Rij) 

Jo Jo 


( 11 ) 


where 7^ x (R,y) is the spatial covariance function for the wet refractivity fluctuations, 
evaluated at 


R ij = r (e,-, 8i, z) - r (tj, 6j, z') + vA t {j 


( 12 ) 


The covariance function can be written in terms of the structure function D x as 


H X (R) = 7l x (0) - -D X (R) 


( 13 ) 


The expression (11), along with the expression (1) for the structure function for Kol- 
mogorov turbulence, enables us to evaluate the T-y and hence the covariance matrix 
of estimated atmospheric parameters from (S). 
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3. Estimating Gradient Parameters: Numerical Results 


We have undertaken to evaluate numerically (8) for the specific equivalent zenith 
delay model 


T (,e,t — T o t- V~(t — t 0 ) 4- (G n cos 6 cot +G e sin 9) cot e (14) 

The model (14) contains: a zenith delay parameter r* representing the zenith delay at 
the reference epoch t 0 ; zenith-delay rate parameter V.; and two zenith-delay gradient 
parameters: G n for the north direction and G t for the east direction. (The azimuth 
angle is measured east from north.) The development of (14) is given in Davis et al. 
[1992]. 

To evaluate (8), we require a “schedule” of observations. We have chosen a sched- 
ule similar to that used by the WVR in the study in Davis et al. [1992]. This schedule, 
shown in Table 1, consists of a 360° azimuth scan at an elevation angle of 30°, followed 
by “tip” curves in the cardinal azimuths. 

In performing the numerical integration, we expressed (11) using (13) and (1) as 

Tu = [ k dz 

J 0 

= l * J W*(0) - \ £ dz f h dz' D(Rjj) (15) 

-*-rt*j:**r ’ 

where we have used the fact that 7Z X (0) is the variance of refractivity fluctuations 
and we have expressed the quantity as cr^ 2 , the variance of the zenith wet delay 

(see TL). 


I 


dz' 7£ X (R ij) 
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In performing the required integration numerically, we must specify the values of 
five quantities: (1) the zenith delay variance ; (2) the strength parameter C; (3) the 
effective height of the troposphere h; (4) the wind speed v\ and (5) the azimuth of the 
wind direction 9 V . For this study, we first adopted “standard” values for , C, and h 
(see Table 2), and performed the calculations for a range of wind speeds and azimuths. 
In the next section, we discuss the sensitivity of our results to changes those values. 

The results of the calculations with the standard values are presented in Figure 1. 
This plot shows the standard deviation, i.e., the square root of the diagonal element of 
(<5p<5p T ) from (8), for the north gradient, G n in (14), as a function of wind direction, 
for a range of wind velocities. The interpretation of this standard deviation is that, 
under the assumptions of Kolmogorov turbulence and the frozen flow model, and 
observations in directions and at times as shown by Table 1, and given the gradient 
model (14), if we have an ensemble of estimates of G n , the standard deviations of those 
estimates will be as shown in Figure 1. (The expectation of the gradient estimates is 
zero.) 

For the case v = 0 m s -1 Figure 1 indicates that there is no dependence on wind 
direction. For higher wind velocities, two major changes occur: the standard deviation 
of the gradient parameter decreases, and the variation with wind direction becomes 
more pronounced, with larger gradient estimates being more probable in directions 
orthogonal to the wind direction. This last result might seem counterintuitive, since 
one normally thinks of the wind as bringing in new air masses which may be more or 
less moist than the air mass currently over the site (as with the passage of a front). 


S 


Our model, however, did not include any real changes in atmospheric conditions, and 
Figure 1 represents the results for an atmosphere whose spatial variations are due 
solely to Kolmogorov-type (zero mean) turbulence. The explanation for the shape of 
the curves in Figure 1 lies instead with the effect of the wind velocity on the R tJ of 
(15). Referring to Figure 1, for a north-south wind, the R tJ - oriented in an north-south 
direction are “longer,” through (10), than those in the east-west direction, for equal 
elevation angles. In other words, the wind tends to effectively “compress” the field 
of fluctuations in the direction of the wind. The result of this compression is that 
the refractive index variations along the direction of the wind are averaged over a 
greater distance, with a greater probability of being averaged out. In this manner, a 
prevailing wind with no true static gradients may result in larger gradient estimates 
for the horizontal direction perpendicular to the wind than for that parallel to the 
wind, and may give the impression of a prevailing gradient. 

This explanation also enables us to understand the asymmetries in Figure 1. 
This plot shows that the standard deviation for G n for a wind blowing to the north 
is much smaller than that for a wind blowing to the south. This results from the 
particular observation schedule we chose, which starts with an observation in the north, 
and continues with observations at increasing azimuth angle. Thus, a wind directed 
towards the north blows the turbulence field away from the immediately subsequent 
observations, whereas a wind directed towards the south blows the turbulence field 
towards the immediately subsequent observations. 
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4. Sensitivity to Variations in Model Parameters 


In the previous section, we presented our “standard” values for various parameters 
used in the numerical evaluation of (SpSp T ). It is straightforward to show that for the 
symmetric observing schedule of Table 1, the variance of the gradient estimates does 
not depend on our choice for a", . Figure 2 shows three curves. The solid line is our 
calculation using the standard values and for v = 8 m s -1 . The other curves represent 
values of C 2 and h increased by 50%. We can see that the size of the standard deviation 
of the estimate of G n is sensitive to the values for these two quantities, but that our 
conclusions regarding the dependence on the wind azimuth are not. In both cases, 
the size of the gradients increase, for the same reason. In both cases the changes to 
the parameters cause an increase in the variance of the delay fluctuations, and larger 
gradient estimates can therefore be expected. 

Figure 3 shows the results of our calculations for a change in the observing sched- 
ule. The elevation angle for the azimuth scan has been changed from 30° to 10°. 
In this case the gradient estimates will be much smaller, since, as with higher wind 
speed, the refractive index fluctuations will be averaged over much larger distances. 
For this wind speed, the direction of the wind has a much smaller effect since the large 
separation between the ray paths causes the correlation between observations to be 
small. 
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5. Ramifications for Microwave Radiometric and Space Geodetic Measurements 


In the preceding sections, we calculated the effects on gradient parameters esti- 
mated from ground-based radiometric data of spatial fluctuations in the wet refractiv- 
ity. We showed that a prevailing wind can lead to larger gradient estimates in some 
directions than in others, even when no prevailing gradients are present. The question 
arises as to whether our method for the analysis of the data was in some way faulty. 
In fact, our model for the equivalent zenith delay observations and our least squares 
inversion method did not take into account either of our basic assumptions concerning 
the structure of atmospheric turbulence and frozen flow. If these assumptions could 
be properly built into the analysis, no azimuthal variations should be apparent in 
Figure 1. These variations are thus understood to result from the limitations in our 
analysis model causing a coupling between the spatial/temporal correlations in the 
refractivity fluctuation field and the observing schedule. 

For example, suppose that our data type were Global Positioning System (GPS) 
phase observations. GPS receivers are able to observe several satellites from different 
elevations and azimuths simultaneously. If we were to estimate atmospheric delay gra- 
dient parameters from each epoch of GPS observations, then it would be similar to 
our results in Table 1 for v = 0 m s -1 . Thus, there would be no preferred gradient 
direction for this single-epoch estimate. However, the correlation among atmospheric 
fluctuations effectively induces temporal and spatial correlations in the GPS observa- 
tions, a case similar to that investigated for very long baseline interferometry (VLB I) 
by TL. 
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In the case of VLB I, observations are obtained at low elevation angles, often below 
10°. In this case, the relevant curve is the e = 10° in Figure 3. This shows that VLB I 
observations will also be less affected in the turbulence-only case by wind, and the 
estimated gradients will probably be much smaller than those estimated either from 
GPS or ground-based radiometric data. Rogers et al. [1992], who compare gradient 
parameters estimated from VLBI and radiometric data, seem to confirm this result, 
however the comparisons are made for different sites. 

V'e have presented here results for an atmosphere which is unchanging except 
for (zero mean) fluctuations described by Kolmogorov turbulence. Clearly, there are 
instances where data clearly indicate the presence of a persistent azimuthal asymmetry 
[e.g., Dixon and Kornreich Wolf, 1990; Davis et al., 1992]. Such asymmetries may be 
caused by horizontal variations in humidity or temperature. The results in this paper 
describe expectations for estimated gradient parameters in the “background” case 
where no such horizontal variations exist. 

Acknowledgements. I would like to thank J. Mitrovica for useful comments. This 
work was supported by NSF grant EAR-9105502, USGS grant 1434-92-G-2170, NASA 
grant NAG5-1930, and the Smithsonian Institution. 
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Table 1. Observation “schedule.” 


Azimuth 

(deg) 

Elevation 

( de g) 

Time 

(sec) 

0 

30 

0 

20 

30 

20 

340 

30 

340 

0 

20 

360 

0 

30 

370 

- 

90 

380 

180 

30 

390 

180 

20 

400 

90 

20 

420 

90 

30 

430 

- 

90 

440 

270 

30 

450 

270 

20 

460 
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Table 2. Standard values for parameters. 


Parameter^ 

Value 


5 mm 

C 

2.4 x 10" 7 m" 1 / 3 

h 

1 km 


! Standard values for C and h taken from Treuhaft and Lanyi [1987]. 
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Figure Captions 


Fig. 1. Standard deviation of the estimated north gradient parameter G n as a function 
of wind direction, for various wind speeds. The observing schedule of Table 1 and the 
standard values for a^ z , C, and h from Table 2 were used, (a) v = 0 m s -1 . (b) 
v = 2 m s -1 . (c) v = 4 m s -1 . (d) v = 6 m s -1 . (e) v = 8 m s -1 . (f) v — 10 m s _1 . 

Fig. 2. Standard deviation of the estimated north gradient parameter G n as a function 
of wind direction, for a wind speed of 8 m s -1 . The observing schedule of Table 1 was 
used, (a) Standard values from Table 2. (b) C 2 = 1.5x standard value, (c) h = 1.5x 
standard value. 

Fig. 3. Standard deviation of the estimated north gradient parameter G n as a function 
of wind direction, for a wind speed of 8 m s _1 . The standard values from Table 2 were 
used, (a) Observing schedule from Table 1. (b) Observing schedule from Table 1, 
except the azimuth scan was performed at an elevation angle of 10°. 
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